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Abstract 

Further investigations of implicit solutions to non-linear partial 
differential equations are pursued. Of particular interest are the equa- 
tions which are Lorentz invariant. The question of which differential 
equations of second order for a single unknown (j) are solved by the 
imposition of an inhomogeneous quadratic relationship among the in- 
dependent variables, whose coefficients are functions of cj) is discussed, 
and it is shown that if the discriminant of the quadratic vanishes, 
then an implicit solution of the so-called Universal Field Equation is 
obtained. The relation to the general solution is discussed. 

Implicit Solutions of PDE's 
1 Introduction 

The study of nonlinear equations particularly those incorporating Lorentz 
invariance has demonstrated that a characteristic feature of large classes of 
solutions, and sometimes the generic solution is that solutions are obtainable 
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only in implicit form. As such, a few implicit solutions to nonlinear equations 
are recorded in the literature, PP [21 jS] 111 [Zj , one of the neatest being 
the result that the general solution of the Bateman equation, 

{4>xf<Ptt - '^<Px4>t4>tx + {.4>tf4'xx = 0, 

where subscripts denote derivatives , is given by the implicit solution of the 
following constraint between two arbitrary functions of one variable; 

tF{(P) + xG{4>) = 1, (1) 

where F, G are arbitrary functions of 0. Since this solution of a second 
order partial differential equation depends upon two arbitrary functions of the 
space variables for a hypebolic equation, and the functional form of {F, G) 
an be fitted to the intitial configuration and time deivative, it is the most 
general. This result is well known. One of the simplest ways of understanding 
it is to look for a solution in the implicit form 

t = A{(P,x), 

Differentiation of this equation with respect to {t, x) and substitution of 
the partial derivatives of (p into the Bateman equation shows that this func- 

d'^A((f) x) 

tion will be a solution provided — = 0. But this simply means that 

ox^ 

A{(p, x) is linear in x, i.e. 

t = Ao{(j)) + Ai{<f))x, 

which is tantamount to (P). Later we shall generalise this approach. The 
Bateman equation possesses a remarkable property of covariance, which is 
reflected in the transformation used to obtain it, and also in form of the 
solution (0), namely that if (p{x,t) is a solution, so is any function of (p. In 
this paper, instead of asking for the solution to a given nonlinear PDE, we ask 
for those equations which are solved by an implicit functional relationship 
among the unknown and the independent variables A natural attempt to 
generalize the Bateman result to the linear constraint 

i=n 

'^Fi{(j))xi = constant, 

i=l 

leads to the surprising conclusion that the function (f) defined thereby is a 
'universal' solution to any equation of motion arising from any Lagrangian 
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dependent upon (f) and its first derivatives, and homogeneous of weight one in 
the latter jU] . It so happens that if there are only two independent variables, 
the equation resulting from this prescription is unique, and is the Bateman 
equation. In the case of more than two independent variables, however, the 
linear constraint cannot be designed to fit arbitrary initial conditions and is 
not the most general solution. One of the aims of this paper is to extend this 
study to quadratic constraints of the form 

i=n j=n 

^^M,,(0)x.x, = 1. (2) 
i=i j=i 

Note that both the linear and quadratic constraints remain form-invariant 
under a general linear transformation of the independent variables Xi, hence 
under Lorentz or Euclidean transformations, so viewed as classes of implicit 
solutions, they should apply to Lorentzinvariant equiations. The quadratic 
system has been studied before, by M. Dunajski ^Ujj from a different point of 
view. He proposes to determine the conditions upon (f) such that an equation 
of D'Alembertian type; 

dxi \^dxj) 

admits a solution on the manifold defined by (|21). A particular example of 
this occurs as an exercise in Jeans' famous book on Electromagnetism [TT] 
'Shew that the confocal ellipsoids 



a2 + A + A c2 + A 

can form a system of equipotentials and express the potential as a function 
of A.' If we take © and differentiate it with respect to x^ the result is 

i=nj=n i=n 

E M[^{(t>)x,x,4>p + E 2Mp,((^)x, = 0. (3) 
i=i j=i j=i 

Here (pp denotes differentiation of (p with respect to Xp, and a prime means 
differentiation with respect to the argument. Differentiation again gives 

M-j{(l))XiXj(f)p(l)q + M-j{(f))XiXj(l)pg 

+2 (m;.((^)x,(^p + MV(</>)x,-0, + Mp,i<j))) = 0, (4) 
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where repeated indices imply summation. Multiplication by r]pq, summing 
and substituting from for (pp etc and integrating gives the equation (2.9) 
of Dunajski's paper 

where 2g' = r]ijM'^K This equation is the starting point to his analysis. 

2 Alternative interpretation. 

Here the purpose is to ask the different question 'What differential equations 
possess solutions of implicit form given by the quadratic ansatz (j2))?' This is 
an obvious extension of the linear situation. Taking equations Q for indices 
{pp,pq,qq) respectively it is possible to eliminate the second derivatives of 
Mij and at the same time many of the other terms. If we define A by 

i=n j=n 

i=i j=i 

and fpq by 

fpq = {<f)pY<Pqq " 2(j)p(J)g(j)pg + {(pgfcppp (S) 

where subscripts denote partial differentiation with respect to Xp etc, so that 
fpq = 0, (jHI) is simply the Bateman equation in the independent variables 
), the resultant eliminant is just 

Xfpg + (Mpp{(t>gf - 2Mpg(j)p(j)g + MggitPpf) = C (Q) 

In fact, there is a more general eliminant, namely 

{X(ppg + Mpg)(j)r(j)s - {Hrg + Mrg)<pp(ps 
- {X(f)ps + Mps)(f)r(t)q + {\(prs + M,,)0p0, = (7) 

Rewriting Q in the form 

\(t>p + MppXp + ^ MpjXj = (8) 

The off diagonal elements of the matrix may then be eliminated from (jHJ 
by using (jH)). We see that (jH)) and © constitute |n(n + 1) linear equations 
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for the elements of the symmetric matrix M. Imposition of one further 
relationship upon the matrix elements, which may be taken as homogeneous 
and nonlinear, will leave a nonlinear equation of second order for 0. In 
the next section this procedure is carried out in detail for two independent 
variables. 



3 Two variable case 

In the case of two independent variables the relevant equations are 

\(t)i + Miixi + M12X2 = 0, 

X(j)2 + M22X2 + M12X1 = 0, 

A/i2 + (Mn(</)2)'-2Mi20i02 + M22(0i)') = 0. 

The easiest way to manage this is to solve these equations for 
Mil, M12, M22, giving 

M - A(Xi0f + X20f02 + Xlfi2) _ .2, ,2 X 2. N 

- + 2X1X20102 + - M01-AX,/12), 



A(xf /12 + Xi0i0i + 
{xj(pl + 2XiX20102 + 3^202 
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A(xi0?02 + a;20i02 - 2:1X2/1: 



M12 = - ") air" o "'^n I "212^ = A^(0i02 + AxiX2/i2)- 
(xf(/)f + 2x1X20102 + a;^02) 

Example 1. The imposition of the homogeneous relationship 

M11M22 - MI2 = 0, 
yields the differential equation 

/12 = 0. 

That this is just the Bateman equation is to be expected, since the re- 
lationship imposed makes the initial quadratic factorize into the square 
of (HJ. 

Example 2. Imposition of the constraint M12 = which unlike the 
previous constraint breaks the Lorentz invariance, unless (xi, X2) are 
considered as light cone co-ordinates yields the equation 

0102X1 + 0102X2 - X1X2/12 = 0. (9) 
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Now the constraint equation Q becomes 



Mii(0)x^ + M22(0)x: 



which gives the general solution for the Bateman equation expressed in 
the variables {xf, X2), which becomes the equation Q when written 
in the original variables {xi, X2). 

A similar result can be found in the three variable case; If the discriminant 
of the quadratic form is set to zero, i.e. 



det 



Mi3 



M12 
M22 
M23 



Mi3 

M23 
M33 



(10) 



which remains compatible with the form-invariance of the ansatz under linear 
co-ordinate transformations, then the resulting equation is what has been 
called by its authors 'the universal field equation', because it arises from 
an infinite number of inequivalent Lagrangians, 



det 







n 
hi 



•■P2 
^12 

^23 



•■P2. 
^13 
^23 
^33 



0. 



(11) 



This result is a special case, shortly to be generalized. 



3.1 Multivariable case 

The simplicity of the solution suggests that it might be possible to solve 
equations (0) and (jHl) in general. Indeed this is so, and the solution is, using 
the result responsible for the simplification above, namely 

Y.X(t)jXj = -1. 

i 

The result is 
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where 



In particular; 



[12) 



The imposition of any relation among the Mpg replacing, if necessary, the 
parameter j hj —J2 4>i^i will provide a second order differential equation 
whose solution is given in implicit form by the original quadratic relation 
(121) In particular, the imposition of the relation det \ Mjk\ = imposes the 
condition that the implicit solution of the quadratic relation (0) for solves 
the Universal Field Equation 



det 







0. 



(13) 



Proof 



Consider the equation 



det 



62 

n 









0? 


01 . 


■ 4>l 





fl2 ■ 


■ fin 


/l2 


. 


■ f2n 


fin 


f2n ■ 


. 



(14) 



From the definition of /p^, (0), by elementary row and column operations 
this equation is equivalent to 



:-l)"2"-^n'/''det 



i=l 







0, 



i.e. the Universal Field Equation. On the other hand, from (jH)) the left hand 
side of equation (fTH) is equivalent to 



r) \ n—l i=n 



0i 
0, Mij 



(15) 



Now if det \Mjk\ = 0, then for a consistent non trivial solution of the other 
n linear equations which determine full set of the matrix elements Mpg in 
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addition to (jH)), namely (jS)), we see that every determinant obtained from 
det \ Mjk\ by replacing each each column in turn by a column vector of the 
must vanish, hence the entire determinant (fT^ vanishes, and the result 
is proved. 



4 Generalisation 

It it has already been observed that the Universal Field Equation is solved by 
a linear assumption [2] [HI, and now it has been shown that a quadratic ansatz 
also works. In addition, explicit solutions can be obtained by taking as a 
homogeneous function of its arguments of weight zero [2]. characterization 
of the general solution, which is known to be integrable, [S] follows similar 
lines to those of 0. 

The clue as to how to do this comes from the discussion in the introduction 
in which solutions of the form 



t — ^4(0, Xj), Xn — t, 



are sought. If this satisfies the Universal field equation, then the homoge- 
neous Monge- Ampere equation must be satisfied. 



det 



dxidxj 



0. 



Both the linear and quadratic relations imposed here satisfy this constraint. 
As a second example, consider the equation arising from the Lagrangian 

+ + (Pl, namely 

/l2 + /23 + /l3 = 0. 

Setting xs = t, this requires as a condition upon A 

(1 + Al)A22 + (1 + ^2)^11 - 2A1A2AU , : 



0, 



the so-called two dimensional Bateman equation, solution methods for which 
are described in PSji^B]- Thus, as in the previous case, the general solution 
depends upon the implicit solution of a similar non-linear equation with one 
less independent variable. 
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5 Complex Bateman equation 



This final chapter is a bit of a diversion from the main theme of the article, 
but is included as it is both anelegant result, obtained here in a slightly 
different way, and a nice example of how results in mathematical physics 
sometimes have been already anticipated several years earlier. The Bateman 
equation possesses a complex form 

(Px(pz4'yw - 4>x4>w4>yz + 4'y4>z4>xw " 4>y4'w4'xz = 0, (16) 

whose general solution has been published in 1935 by Chaundy [T] , redis- 
covered and generalized in 1999. P3I. It is given by equating an arbitrary 
function of {x, y, (f)) to another of {z, w, (p), i.e. 

Fix,y,(f)) = G{z,w,(p). 

This equation may be expressed in first order form; 

dv dv 
dz dw' 

— = V— (19) 

— = u— (20) 
dz dw' 

The way this works is that we have 

F F 



G(f, — F^ Gfj, — F^ 

Gin 



Hence, from (jl9p and (j2(Jj) we have 

Fx Gz 

v = y; u = — 

J. y 

i.e. 

5m _ / Gj^ _ Gjn^Gj\ ^ _ ^du _ ^du 
dx \G^ Gl ) <Pydy dy' 
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thus verifying (fT7|) . and similarly for (fTS|) . 

There is a huge class of explicit solutions to the equation (fTHj) . Take any 
two arbitrary differentiable functions, f{x,y) and g{z,w) and construct an 
arbitrary function F{f{x,y),g{z,w)). Then 

(f) = F{f{x,y),g{z,w)) 

is a solution to (fTB|). 
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